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Spectroscopic amplitudes play an important role in nuclear capture reactions. 
These amplitudes are shown to include both single-particle and polarization effects: 
the former through their spatial dependence and the latter through their normaliza- 
tion (the spectroscopic factors). Coupled-channels equations are developed for the 
spectroscopic amplitudes. These equations serve as a convenient starting point for 
the derivation of several approximations: Hartree, Hartree-Fock and two different 
single-particle models. The single-particle models include antisymmetry in different 
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ways, but both miss many-body effects. Therefore, cross sections calculated with 
either of these models need to be multiplied by the spectroscopic factor. 

PACS numbers: 25.40.Lw, 24.50.+g, 21.60.-n, 26.65,+t 



I. INTRODUCTION 

X 

Nucleon capture reactions at low energies, such as 7 Be(p, 7) 8 B, 16 0(p, 7) 17 F*, or 
7 Li(n,7) 8 Li, play an important role in our understanding of astrophysical phenomena. For 
example, in the hydrogen-burning process in stars such as our sun, low-energy proton cap- 
ture on beryllium takes place in the proton-proton chain, and the 16 0(p, 7) 17 F* reaction 
occurs in the CNO cycle [p], §, [§. Exact knowledge of the reaction rates is necessary for 
modeling the energy generation and evolution of hydrogen-burning stars. In addition, the 
7 Be(p, 7) 8 B reaction at solar energies (E cm < 20 keV) plays a key role in the 'solar neutrino 
puzzle' since the neutrino event rate in the existing chlorine and water Cerenkov detectors is 
dominated by the high-energy neutrinos produced in the subsequent (3 decay of 8 B [I], @, f| . 
The 7 Li(n, 7) 8 Li reaction is a key element of primordial nucleosynthesis in inhomogeneous 
big bang scenarios jT|, H . It initiates a sequence of reactions which bridge the mass gap at 
A=8 and thus its rate is crucial for determining the amounts of heavier elements produced 
in these models. 
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Direct measurements of capture reactions at energies corresponding to astrophysically 
relevant temperatures, however, are often very difficult, since the cross sections diminish 
exponentially at low energies. Thus, theoretical studies of these processes become very 
valuable. In addition to microscopic theories, such as the nuclear shell model or cluster 
models, one-body potential models provide a popular framework for such investigations. 
For example, the potential model was used in ref. |J to discuss the energy dependence of 
the reaction rates. In the one-body potential model, a single-particle wave function is used 
to calculate various observables; microscopic substructure effects are partially accounted 
for through the use of spectroscopic factors. This strategy, however, has recently been 
the subject of vivid discussions [^ ||. At issue is the proper normalization of the cross 
section. Csoto J7|, for instance, maintains that spectroscopic factors should not be included 
in potential-model calculations of the 7 Be(p, 7) 8 B cross section at low energies, since the 
reaction depends only on the asymptotic normalization coefficient of the 8 B bound-state 
wave function |J, whereas the spectroscopic factor arises from the short-range properties 
of the wave function. Mukhamedzhanov et al. [||, ^ argue in favor of using a different 
approach, based on asymptotic normalization coefficients instead of spectroscopic factors, for 
determining the relevant cross sections. However, the asymptotic normalization coefficient 
contains short-range effects. It can actually be given in terms of an integral over the interior 



of the nucleus 110] and its interdependence with the spectroscopic factor has been noted in 



the earlier work of Locher and Mizutani |IT| and of Lovas et al. fll2 

These recent discussions have motivated us to revisit the question of the proper treatment 
of microscopic nuclear structure effects in one-body models. In the present work we focus on 
the role of spectroscopic amplitudes and factors. The use of spectroscopic factors in nuclear 
reaction calculations dates back to the early days of nucleon transfer reactions |T3f and 



continues to be central in the interpretation of such processes |TJj . With the renewed interest 
in nucleon capture reactions in the context of astrophysical scenarios, and the emerging need 
for very accurate reaction rates, it becomes necessary to review and clarify the assumptions 
associated with one of the most basic models of nuclear physics, the one-body potential 
model. 

Before proceeding with the formalism, it is necessary to clarify the meaning of the terms 
spectroscopic amplitude, spectroscopic factor, and one-body model. For example, it is im- 
portant to realize that there are many spectroscopic amplitudes and spectroscopic factors 
associated with an A-body system, namely one corresponding to each excited state of the 
(A — l)-body system. The different amplitudes are not independent since they are related 
by a model- independent sum rule fL5| . In addition, the spectroscopic amplitudes are related 



by a set of coupled-channels equations. 

The question of what is meant by a one-body model is more complicated since different 
notions are associated with the term. The coupled- channels formalism presented here is 
used to derive two different one-body models and to study the connections between them. 
One approach, which uses either the Hartree or the Hartree-Fock approximation in an in- 
termediate step, leads to A noninteracting nucleons in a one-body potential. This, strictly 
speaking, is not a one-body model since there are still A particles. However, the orbitals 
corresponding to the different particles decouple and the problem reduces to solving A one- 
body equations. The fact that we still have A nucleons explicitly present has the advantage 
that antisymmetry can be built in ab initio using a Slater determinant for the A-body wave 
function. Another method for generating a one-body model is based on integrating out the 
coordinates of (A — 1) nucleons, effectively projecting onto a low dimensional space. In 



this approach antisymmetry is enforced separately for each channel in the coupled-channels 
equations. Truncating to a single channel results in the approximation advocated by Varga 
and Lovas et al. [IB|, who study the cluster substructure of 6 Li in this framework. Their 
work is based on the generator coordinate method. The derivation presented here serves to 
clarify the relation between their model and other approximations. 

In Section II, we define the terms spectroscopic amplitude and spectroscopic factor, derive 
coupled-channels equations for the amplitudes, and discuss several approximation schemes, 
including the first one-body model described above. An alternative approach to the spectro- 
scopic amplitudes is considered in Section |T|, which includes details on the second one-body 
model mentioned above. In Section fV], we explore the physical aspects associated with the 
spatial dependence of the spectroscopic amplitude and its norm, the spectroscopic factor. 
We give an expression for the reaction matrix elements in terms of the spectroscopic am- 
plitudes. For both single-particle models considered here the reactions rates need to be 
multiplied by a spectroscopic factor. Our conclusions are summed up in Section [V|, and 
various technical aspects of our work are included in the Appendices. 

II. SPECTROSCOPIC AMPLITUDES 

In this section we expand an A-body wave function in terms of the ground and excited 
states of the (A-l)-body system and derive a set of coupled differential equations for the 
expansion coefficients, of which the spectroscopic amplitudes are a special case. The resulting 
equations of motion give insight into the long-range behavior of the A-body ground state 
wave function and lead to various approximation schemes. In particular, the Hartree and 
Hartree-Fock approximations and a one-body potential model approach emerge naturally. 
For simplicity, spin and isospin degrees of freedom are suppressed and the Coulomb potential 
is neglected. Furthermore, we do not consider the center-of-mass motion here; this subject 
is discussed in Appendix |B|. 

Let \E f ^_ 1 (ri, • • • , Ta-i) denote the n-th (fully antisymmetric) eigenstate of the (A — 1)- 
body Hamiltonian Ha-i- The collection of all \l^4_i> including both bound and continuum 
states, forms a complete set of states in the space of antisymmetric (A — l)-body wave 
functions. Using this set, {^A-i}n=i,2,..., one can construct a basis for the space of antisym- 
metric A-body wave functions, by defining ^ A ,r (ri, ■ ■ ■ , Va) = A.^' r (r*i, • • • , 7*a), where 
\l/^' r (ri, • • ■ , Va) = ^A-i{. r ii " • • i r A-i)8(r — ta) and A antisymmetrizes between the A-th 
coordinate, which occurs in the delta function, and the (A — 1) coordinates in ^a-v This 
'intercluster' antisymmetrization operator is normalized so as to satisfy A 2 = yAA. 

The \l/^' r span a space that includes both totally antisymmetric A-body states and mixed- 
symmetric states, which are antisymmetric in the first A-l coordinates and symmetric under 
interchanges involving the A-th coordinate. The basis sets {^a-i}^ {^a^}; an d {./i^a^} are 
discussed in more detail in Appendix [A|. 

An arbitrary antisymmetric A-body wave function ipA can then be expanded as: 
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with expansion coefficients: 



0»(r) = U*2 IVU> 

= >/I<*7|V*> (2) 

Here a(r) [a^(r)] is an annihilation [creation] operator which destroys [creates] a nucleon at 
position r and obeys the usual anticommutation relations. Note that only the first expansion 
in eq. ([!]) is manifestly antisymmetric, while in the other two expressions the antisymmetry 
information is carried by the expansion coefficients. The coefficients are identical in all 
three expansions if and only if if) a is fully antisymmetric. In this case, the mixed-symmetry 
components of ^ n A ,r do not contribute. 

When if) a denotes a bound state, the <f> n (r) are called spectroscopic amplitudes, and the 
associated integrals 

S n = fdr\cf> n (r)\ 2 (3) 



are the frequently used spectroscopic factors jr7|. They obey the sum rule Y^Li ^n = 
A (see ref. [[HJ and also eq. ( |A21| )). We observe that for a bound A-body state, there 
are many spectroscopic amplitudes (and thus many spectroscopic factors), namely one for 
each excited state of the (A — l)-body nucleus. Given the structural information on the 
(A — l)-body system that enters the wave functions ^a-u the spectroscopic amplitudes 
completely determine the wave function if) a- The spatial dependence of (f) n { r A) is related to 
the properties of the single-particle orbital of the A-ih nucleon in the larger system, and the 
norm of <f> n {r), the spectroscopic factor, provides a measure of the structural similarity of 
the n-th excited (A — l)-body state and an (A — l)-body subcluster of the full ^4-nucleon 
system. 

The interdependence of the spectroscopic amplitudes can be made explicit upon deriving 
a set of differential equations for the (f> n (r). For an A-body Hamiltonian Ha, which contains 

a kinetic energy term and a two-body potential, Ha = — J2i=i 2^T" + f Si ?=i ^(l r * — r il)' we 

can write H a = Ha-i — ^f- + 12i=i V( I r i ~ r A I ) , where Ha-\ denotes the Hamiltonian of the 
(A — 1) body system and m, is the mass of the i-th nucleon (for simplicity we assume m, = m 
for i — 1, . . . , A). Since |^_i) and \i])a) are eigenstates of Ha-\ and Ha, respectively, we 
have: 
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(n-i\H A \ipA) = E A {n.M 

Inserting the expansion given in eq. ([!]) for \iPa), we obtain a set of exact, Schrodinger-like, 
coupled equations for the spectroscopic amplitudes: 

_2 oo A— 1 

(E A - El_AUr) = ~^Mr) + £<*5-il E V ^ ~ r \)\*A-i)Mr). (5) 

n'=l i=l 



This set of equations, originally derived for stripping reactions (see for example ref. JT8|), is 
not sufficiently appreciated. In the form just given, these equations may be too complex for 
use directly in calculations, since the V nn >(r) = {®a-i\ ^2i=i ^(l r * — r l)l^!4-i) couple an 
infinite set of coefficients. However, they allow us to discuss the long-range behavior of the 
A-body ground state wave function as well as derive commonly used approximations. 

The solutions of eq. (f|) include not only the physically relevant completely antisymmetric 
states, but also the mixed- symmetric, i.e., unphysical, states. As long as the coupled equa- 
tions are solved exactly this does not cause any problems since the two classes of states do 
not mix. However, when approximations are invoked one needs to ensure that the unphysical 
solutions are eliminated. 

From Eq. @ we can extract information on the long-range behavior of the v4-body 
ground state wave function. The matrix element V nn /(rA), for n or n' corresponding to an 
(A — l)-body bound state, has a range which is determined by the convolution of the (short- 
range) two-body potential with the bound-state wave function. Thus it falls off rapidly 
with increasing ta = \r^\. When both ^a-\ anc ^ ^a-i describe continuum states, the 
matrix element has a long-range, but infinitesimal, tail. Solving eq. (|^) outside the range 
of the potential, we find that the 4> n iTA) decouple and fall off exponentially, 0«(^a) oc 
exp[— ^j2m{E\_ x — E A ) r A ]/r. Since (E\_ x — E A ) is smallest for n = 1, the long-range 
behavior of the A-body wave function is dominated by (f>i(rA)'&A-i( r i7 ' ' ' > r A-i)- Thus 
processes with reaction probabilities that are peaked in the asymptotic region depend only on 



the asymptotic normalization coefficient, A nc = lim^oo 0i(r)r expL 2m(E A — E\_^)r] (see 

also refs. || |9]). The asymptotic normalization coefficient is a property of the spectroscopic 
amplitude and hence implicitly includes the spectroscopic factor. 

To obtain an explicit expression for the asymptotic normalization coefficient, we Fourier 
transform eq. ([5]) (to include center-of-mass corrections use eq. ( pil|) ). This gives us: 
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Mk) = ~ W E / dr cttr W-il E V ^ ~ r \)\*A-x)Mr). (6) 

Ea - E A _ X - ^ n , =1 J . =1 

The pole in <ft n (k) is thus explicitly seen. It is this pole which is responsible for the upturn 
in the astrophysical S factor of the 7 Be(p, 7) 8 B and 16 0(p, 7) 17 F* reactions as the incident 



momentum goes to zero 0. As shown in ref. |L0|, the residue at this pole for n = 1 is 



proportional to the asymptotic normalization coefficient. We have: 

A-l 



A nc = \unjm7r 2 E / dr e ife "<^_ilE y (^ " r \)\<-i)Mr) , (7) 
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where k = JjImiEA ~ -^a-i)' i- e - the asymptotic normalization coefficient can be given 

as an integral over the interior of the nucleus (since ^l\-\ nas the spatial extent of the 
ground state of the (A — l)-body system). The price paid is that there is a sum over all the 
spectroscopic amplitudes. 

For arbitrary distances, we may consider approximating the ^a-i i n ec i- © by Slater 
determinants constructed from the spectroscopic amplitudes 0»(r) with i G {1, . . . ,A}. In 
this case, we obtain the Hartree-Fock equations in coordinate space. The A Hartree-Fock 
orbitals can thus be identified as approximations to the spectroscopic amplitudes and the 
Hartree-Fock single-particle energies approximate single-nucleon separation energies. The 



local (Hartree) term arises from the diagonal terms (V nn ) in eq. (|^), and the off-diagonal 
terms (V nn i with n ^ ml) give the nonlocal (exchange) potential, i.e., the Fock terms originate 
from the channel coupling. If the off-diagonal elements (V nn i, n ^ nl) are neglected we obtain 
the Hartree approximation. 

A one-body potential model, which treats the nucleus as a system of A noninteracting nu- 
cleons, can be obtained from either the Hartree or the Hartree-Fock approximation. Starting 
from the Hartree equations, one can express the many-body wave functions ^a-i as Slater 
determinants of the <f)i(r), i 6 {1, . . . , A}, and impose the condition that the one-body po- 
tential be the same for each single-particle orbital, V nn (r) =>■ V(r). Alternatively, beginning 
from the Hartree-Fock equations, one can dictate the replacement: 

A-l 

V nn ,{r) = (rk-il^Vdn-r^^i^^Vir)^ . (8) 

In both cases, one obtains the following equation for the (j) n (r): 

(E A - E n A _Mr) = ~Mr) + V{r)4> n {r) , (9) 

which defines a one-body potential model with A non-interacting nucleons in a common 
potential. Thus, the one-body potential model may be regarded as an approximation to 
either the Hartree or the Hartree-Fock approach. In this model, antisymmetry can be 
ensured in a straightforward manner: many-body wave functions, such as ^a-i or if) a, are 
- like in the Hartree-Fock picture - Slater determinants constructed from A given single- 
particle wave functions. 

The Slater determinants, in one-body potential models, usually play only a formal role 
and disappear from sight in actual calculations, to the extent that their existence is of- 
ten forgotten causing conceptual confusion. For the matrix elements of a one-body op- 
erator, 0(r), between states that differ only in one orbital, the remaining (A — l)-body 
orbitals integrate out leaving just the active orbitals. Thus we have a matrix element like 
J dr<j) m (r)0(r)<f) m t(r) where m and ml denote the active orbitals. At first sight this expres- 
sion is a pure one-body expression that seems to have no antisymmetry present and indeed 
at the computational level this is true. In calculating the matrix elements of one-body op- 
erators there is no need to explicitly consider the Slater determinants. However, the Slater 
determinants did their job of ensuring antisymmetry before they were integrated out. 

Since the single-particle Schrodinger equation, eq. (||), originates from the coupled- 
channels equations for the expansion coefficients, eq. (||), the bound-state single-particle 
wave functions of the potential model discussed here approximate the spectroscopic ampli- 
tudes. Therefore, the potential-model approximation gives A nonzero spectroscopic ampli- 
tudes, all of which are normalized to one. Consequently, the associated spectroscopic factors 
are one and carry no information on the many-body correlations of the nuclei involved. 

One can, however, move beyond the simple picture of A noninteracting particles in a 
common potential and include some many-nucleon correlations in a schematic manner. As 
a first correction, one usually allows for single-particle orbitals which are different for the 
(A — l)-body and A-body systems. As a result, the expansion of if) a, eq. (^), will have 
contributions from more than A terms. Consequently, the spectroscopic factors will be, on 
average, less than one, and most individual factors will be smaller than one. 

Another important correction to be taken into account involves the treatment of the 
center-of-mass motion and leads to the introduction of an intrinsic spectroscopic amplitude 



(see Appendix 0). In both the Hartree-Fock and the one-body potential- model approaches, 
the center of mass is erroneously confined in a potential. In potential models which are 
based on a harmonic oscillator potential the center-of-mass corrections can be taken into 
account exactly. In this case, the spectroscopic amplitudes corresponding to valence shell 
orbitals increase by a factor [A/ (A — l)] m ' 2 , where m denotes the major oscillator shell under 



consideration, while the amplitudes associated with lower shells decrease accordingly |L5 
The sum rule for the spectroscopic factors remains unaffected. The center-of-mass correction 
causes some of the spectroscopic factors to be greater than unity. 

The above considerations illustrate that information on the structure of the A-body wave 
function of a nuclear system, which is lost in a simple one-body potential model approach, 
can be in part recovered through the appropriate correction to the single particle wave 
functions, that is by multiplying by the spectroscopic factors. In Section |V| we will discuss 
how the resulting deviations of the spectroscopic factors from one (or zero) affect physical 
observables such as the cross sections for external capture reactions. 

III. AN ALTERNATIVE PERSPECTIVE 

In this section we approach the spectroscopic amplitudes from a different perspective. 
Motivated by cluster-model results, we derive a set of coupled integral equations. While 
general cluster-model calculations involve two or more clusters of arbitrary size, we restrict 
ourselves to describing an A-body nucleus as an (A-l)-body cluster plus a single nucleon. 
The approach pursued in this section has the advantage that explicit channel coupling is no 
longer required to ensure antisymmetry. Again, various approximation schemes can be used 
to simplify the coupled equations. In particular, in one of these schemes the Hartree-Fock 
equations can be recovered; in another a one-body model emerges, which was previously 



discussed by various authors |16], [19| and which differs from the model presented in the last 
section. 

Starting with the Schrodinger equation for an A-body system, Ha4>a = Ea^a, an d 
expanding i\) A i n terms of the antisymmetric basis states jC&Y> we obtain: 



J2 / dr'{ A ^Y\HA\A< y )Mr') = E A J2 dr'( A ^/\ A ^ r ')Mr') , (10) 

n'=l J n'=l J 

or, equivalently: 
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J2 / dr'^(n,rX,r')<Mr') = E A J^ / cW(n, r,7i',r')<MO , (11) 

n'=l "" n'=\ ^ 

where A H(n,r,n',r') = (^YIHaIa^a'^) and A/"(n, r, n', r') = (a^YU^a'^) denote the 
kernels of integral operators XHnn' and J\f nn >, respectively. The operator l nn i (I nn > = j(Hnn> 
or J\ nn '), which is the (n,n') entry of an infinite-dimensional matrix I (I = A 7i or AT), 
acts as follows: I nn >f(r) = J dr'I(n,r,n',r')f(r'). The matrix Jx serves as the norm 
operator for the antisymmetrized basis { A ^Y } and has many interesting properties which 
are discussed in Appendix [FJ. Jx/A is a projection operator which projects an arbitrary 
many-body state into the set of completely antisymmetric wave functions. Thus, eqs. ( p!U|) 
and ( PI ) contain an explicit projection onto completely antisymmetric states, in contrast to 



the coupled differential equations discussed in the previous section, eq. (|5|). This projection 
means that the solutions of eq. (|TTD, (f>n° (r), can contain arbitrary contributions from 
mixed-symmetric states while still satisfying the equations of motion. Thus the projected 
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r'), rather than 0„ (r) corresponds to the 



state, <p n = (1/A) J2n'=i I dr'N(n, r, n', 
physical state. 

With eq. ([[(]) (or, equivalently, eq. (|TTD ) , we have derived a set of coupled-channels (in- 
tegral) equations, analogous to the system of coupled differential equations in the previous 
section, eq. (0). In the present approach, however, antisymmetry is explicitly enforced, 
resulting in a set of equations which contain a projection operator and a complicated ex- 
pression for the effective Hamiltonian. To illustrate this, we compare the Hamiltonian kernel 
in the nonantisymmetrized basis, Tt(n,r,n',r'), with the corresponding expression in the 
antisymmetrized basis, XH(n,r,n',r'). The former is local, i.e., diagonal, in the spatial 
coordinate, although not in the discrete variable n, 



H(n,r,n'y) = ^ n /\H A \^Y) 



5{r 



<L 



h A-l 



2m 



)+(n-i\ y £v(\r i -r\)\ri- 1 ) 



(12) 
(13) 



whereas the latter has off-diagonal contributions from both r and n: 
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A H(n,r,n',r') = J^ dr" H(n, r , n" , r")N (n" , r" , n' , r') 
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^i_ 1 (r,r 2 ,...,r A . 



(15) 



Imposing antisymmetry via the projection operator associated with J\f(n,r,n',r') leads to 
exchange terms in the effective potential. 

For the special case of a one-body Hamiltonian, H = X]i=i-^( r *)> we obtain the 
Hamiltonian kernels H(n,r,n',r') = 5 nn >5(r — r') (E A _ l + H(r)) and /H(n,r,n',r') = 
J\f(n,r,n',r') \E\_ l + H(r)), and eq. ( |i"l"l) reduces to the single-particle Schrodinger equa- 
tion. The kernel J\f(n,r,n',r') guarantees that only expansion coefficients originating from 
an antisymmetric A-body wave function are considered, i.e. M(n, r, n', r') ensures that the 
associated single-particle orbitals are not among the occupied states of the (A-l)-body sys- 
tem. This can be trivially taken into account in the calculations. Nevertheless, even in a 
case as simple as this one, Af(n, r, n', r') is not diagonal in n, n' . 

In order to facilitate working with the above (exact) set of coupled integral equations, 
eq. (0), various approximations may be considered. For example, the Hartree-Fock equa- 
tions are recovered by taking the (A — l)-body basis states, ^_i, to be Slater determinants 
constructed from the expansion coefficients <p n (r). In contrast to the approach presented 
in the previous section, the current method does not require channel coupling to obtain 



the Fock terms. Instead, the nonlocal (exchange) terms are now explicitly present in the 
effective potential, as can be seen in the last line of eq. (|15]). Thus, the advantage of using 
antisymmetric basis states - - coupled channels are not needed to include the Fock term 
contributions --is offset by additional complications in the resulting equations of motion: 
the Hamiltonian kernel is no longer local and the antisymmetry operator appears explicitly 
in eq. (11). Note also that in the Hartree-Fock approximation the channels are implicitly 
coupled through the use of expansion coefficients in the (A — l)-body Slater determinants. 
Another approximation method leads to equations which were previously obtained by 
Varga and Lovas et al. [16|] in the framework of the cluster model. In this approach, we 
ignore those terms in XH(n,r,n',r') and Af(n,r,n',r') that couple different values of the 
discrete variable, i.e. contributions for which n ^ n' holds. The equations of motion for the 
coefficients (f)n{ r ) then take the following form: 

j dr' A H(n, r, n, r')<f> n {r') = E f dr'M{n, r, n, r')0„(r) , (16) 

that is, we have effectively integrated out the coordinates of (A — 1) particles to obtain 
a set of one-body equations. This equation is more general than Hartree-Fock, but the 
Hartee-Fock approximation can be recovered by taking the (A — l)-body wave functions to 
be Slater determinants composed from the spectroscopic amplitudes. We can then derive the 
one-body potential model, eq. (^), as in the previous section. In these cases antisymmetry 
can be exactly imposed through the use of Slater determinants and the explicit projection 
operator, JV, is not needed. 

A second one-body equation can be obtained by taking the (A — l)-body wave functions 
from an external source independent of the spectroscopic amplitudes. This can go beyond 
Slater determinants and Hartree-Fock since the (A — l)-body wave functions can include 
many-body correlations, in principle they can even be exact. The quality of this approxima- 
tion will depend on how well the (A — l)-body wave functions are chosen. In this approach 
the orbitals can be considered one at a time. The price paid for this convenience is twofold. 
First, the expressions for the Hamiltonian become more complicated, compare eq. (|T3|) and 
eq. ([0|). Secondly, when we move beyond Slater determinants, N is explicitly required and, 
moreover, its diagonal element, J\ nn , is no longer a projection operator (see Appendix B), 
so antisymmetry is not explicit. 

The one-body models associated with eqs. (0) and flTEp differ, but in most situations 
both lead to a decoupling of the expansion coefficients. The exceptions are the Hartree-Fock 
approximation and other self-consistent models in which the (A — l)-body wave functions 
are constructed from the expansion coefficients. While the decoupling is in many ways 
advantageous, it also implies that we have lost the information on the relative normalizations 
of the (f) n {r). If we assume that A of the coefficients are nonzero we may take them to be 
individually normalized to one. This is consistent with the overall normalization condition 

In the work of Lovas et al, eq. (|16D was considered for the case n = n' = 1. The Hamil- 
tonian operator h = (An) -1 / 2 aH-u (An) -1 / 2 and the wave function xi — (Aii) 1 ^ 2 ^ were 
defined. These redefinitions have the advantage of leading to an equation which takes the 
same form as the standard Schrodinger equation, h\i = E\i- The \i provides a better ap- 
proximation to the true spectroscopic amplitude than the model 0i (the solution of eq. flTE|)) 
since it includes the effect of (A/n) 1 / 2 , an approximate projection operator onto antisymmet- 
ric states. As noted previously, a projection operator may be needed and this is amplified 
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in the discussion of ref. |7[] in the next section. It is also xi (normalized to one) that in this 
approach should be identified with the single-particle wave function, both because of this 
better correspondence with the true spectroscopic amplitude and because of the form of the 
equation which it satisfies. 

IV. SPECTROSCOPIC AMPLITUDES AND REACTION RATES 

In this section we discuss how the spectroscopic amplitudes can be used to calculate 
reaction rates and elaborate on their physics content. We show that the spatial dependence 
of the spectroscopic amplitude and its norm, the spectroscopic factor, describe different 
physical aspects of the many-nucleon system, namely single-particle properties of the A-th 
particle and distortions of the {A — l)-body system, respectively. We make the connection to 
the potential-model approach and find that both aspects are included, at least approximately, 
when spectroscopic factors are employed to scale the wave functions. Since our findings 
contradict the conclusions of ref. f7j, we explore the claims made in that paper. 

To calculate the relevant reaction cross sections, we expand both the A-body bound 
state, |Vu(**i> • • • j v a)), and the wave function in the incident channel, \ip^ (t*i, •• • , r A )), 
as in eq. flU). Here K specifies the asymptotic momentum of the incident particle relative 
to the (A-l)-body target nucleus. The corresponding expansion coefficients are given by 

<f> n (r) = VA(*Y\Mri, * ' ' > r A )) and <f$(r) = v^ZKK • • • , r A )), respectively. The 
expansion coefficient 4>\{v), associated with the continuum wave function, is the optical 
model wave function |20|, [H]]. Thus the formalism based on expansion coefficients is suffi- 
ciently general to include both spectroscopic amplitudes and optical model wave functions. 
It will be useful whenever we are dealing with one-body operators. The matrix element for 
the one-body transition operator 0(r) can then be written as: 

A 

M = (1>A\520(ri)\1$) = A(1> A \0(r A )\1$) (17) 

= Y.J dr<p n (r)0(r)<j£(r). (18) 



n=l 



To make the connection with the potential model we separate the spatial dependence of the 
spectroscopic amplitude and its normalization, a/S^, as follows: 



S n <Pn(r) , (19) 

where j dr\(j) n (r)\ 2 = 1. In Section II, we have shown that the potential-model wave func- 
tions approximate the spectroscopic amplitudes. Since the norm of <j> n (r) is one, whereas that 
of <p n (r) is S n , the n (r), rather than the <fi n (r) should be identified with the potential-model 
wave functions. There is no equivalent normalization factor associated with the scattering 
state; if)% is normalized asymptotically. 

The transition matrix element given in eq. (|18"D can now be written as M. = 
Er=iV^/dr0;(r)C>(r)0f(r). Both the 7 Be(p,7) 8 B and 16 0(p,7) 17 F* reactions at 
threshold are peripheral, i.e. the capture processes take place at large distances from the 
center of the target nucleus, which is in its ground state. In such situations, as for all direct 
capture reactions, only the first expansion coefficient for ip% contributes and the matrix 
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element reduces to M. — > y/S~l J dr<pl(r)0(r)(j)f(r). Since the cross section is proportional 
to (A^ | 2 , the spectroscopic factor associated with the ground state of the (A-l)-body system, 
Si, occurs linearly in the expression for the reaction rate. Thus we conclude that when a 
potential-model wave function (or any other function normalized to one) is used to calculate 
the transition matrix element Ai, the resulting cross section needs to be multiplied by the 
associated spectroscopic factor in order to account for (some of) the many-body correla- 
tions in the nuclei involved. The implications of this for the asymptotic normalization are 
discussed at the end of this section. 

The separation introduced in eq. (|I9D is motivated by the realization that the spatial 
dependence of the spectroscopic amplitude </> n (r\4), and its norm, the spectroscopic factor 
S n , describe different physical properties of the nuclear many-body system. The former is 
related to the shape of the single-particle orbital of the A-th nucleon in the system, and 
can therefore be expressed through the normalized amplitude (f> n (r). The latter provides a 
measure of the structural similarity of the n-th excited {A— l)-body state and an (A— l)-body 
subcluster of the larger system. Equivalently, the set of spectroscopic factors associated with 
the expansion of ip A , eq. (|l|), can be viewed as describing the distortion of the (A — l)-body 
core due to the presence of an extra nucleon. This can be seen, for example, by casting S n 
into the following form: 

S n = J dr\ct> n {r)\ 2 (20) 

'A-\ \ „ M-l \ 

A 'J I iTri* ( m „ . W, 



AjiTldrAj\Udri)nU(ru---,r A ^W J 



x 



dri/j* A (r u • • ■ , r A _ l} r)4> A (r[, ■ ■ • , r^_ l5 r) 



(21) 



In the last line (expression in square brackets) we have integrated out the dependence on 
the A-th particle. We are left with expressions involving wave functions of the (A — l)-body 
system. The extra particle's influence is still present in the modification it has induced 
in the (A — l)-body cluster. Upon decoupling the equations of motion for the expansion 
coefficients, as required in our derivation of the one-body models, this information on the 
last particle's influence is lost, the many-body correlations contained in the integrals S n 
disappear and the spectroscopic factors become one. 

The spectroscopic factors carry information both on the dynamical distortions induced by 
the interaction between the (A-l)-body system and the extra particle and on antisymmetry 
effects. Specifically, as shown in Appendix y, the spectroscopic factor can be written as: 

S n = N n {^\^Y) , (22) 

where (V'aIV'a ) reflects the influence of the distortions, and the normalization factor iV„ 
contains antisymmetry effects. Specifically, N n keeps track of the requirement that the 
nucleon orbitals in the n-th excited (A — l)-body state must be orthogonal to the orbital 
n (r) of the Ath particle (for details see eq. flU5p). When center-of-mass corrections are 
taken into account, N n may be greater than one, otherwise it is less than or equal to one. 
The matrix element (WaI^a) ^ s always less than or equal to one. Consequently, when 
center-of-mass corrections are ignored, we have the restriction S n < 1 for the spectroscopic 
factor. The influence of these corrections are discussed in detail in Appendix 0. 
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At this point we would like to reassert that it is correct and necessary to include spectro- 
scopic factors in potential-model cross sections of nuclear capture reactions. This procedure 
was questioned in ref . [7| . The author of that paper argues that the short-range correlations 
contained in the spectroscopic factor should modify the bound-state wave function only in 
the nuclear interior and have no effect on the asymptotic behavior. Since multiplying the 
potential model wave function by \fS^, however, affects its overall normalization, including 
in the tail region, the usual procedure of treating microscopic correlations in the potential 
model, viz. through spectroscopic factors, is pronounced to be incorrect. To illustrate his 
point, the author compares a wave function Xc( r )i which describes the relative motion of 
7 Be and p in 8 B, to the spectroscopic amplitude function of the 7 Be+p configuration in 8 B 
[Note that this \ c is different from the wave function xi introduced by Lovas et al. fl6 



and discussed in the previous section]. The wave functions \ c considered individually do 
not contain Pauli effects but, when used in cluster-model calculations, appear behind an 
antisymmetrization operator. The spectroscopic amplitudes, on the other hand, are cal- 
culated from properly antisymmetrized wave functions. The two functions are shown to 
agree with each other and with the appropriately normalized Coulomb- Whittaker function 
in the asymptotic region, but they differ at small radii, as can be seen in Fig. 2 of ref. J7|. 
Their difference is interpreted as a measure of the Pauli effects. Since multiplying Xc{ r ) 
by a spectroscopic factor would affect both its short-range and asymptotic behavior, the 
author concludes that this cannot be the proper procedure for incorporating microscopic 
substructure effects in potential model calculations. 

The argument presented in ref. J7J is not correct. The function Xc(f) does not, in general, 
correspond to a potential- model wave function. First, unlike the one-body wave functions 
we have considered, it already includes the spectroscopic factor through the normalization 
of the A-body wave function. Moreover, very little physics can be associated with \ c outside 
the context of its usual use as a component in a properly antisymmetrized cluster-model 
wave function. To show this, we consider an arbitrary product state ipJKri, ■ ■ ■ ,Ta) — 
(p(rA)ipA-i{ r ii ' ' ' i r A-i)- In the nonantisymmmetrized basis, it has expansion coefficients: 

#(r) = y/A{9YM) = <P(r)VA(n-M-i) , (23) 

and the associated antisymmetrized, i.e. physical state, has expansion coefficients (f) n (r) = 
(l/A) J2ri=i f drj\f(n, r,n', r')0^,(r). In the product given above, ip(r) can be taken to 
correspond to the relative-motion function Xc( r ) of ref. 0. As discussed in the previous 
section, A/" enforces antisymmetry by projecting onto a completely antisymmetric state. 
Thus, the choice of <p(r) is somewhat arbitrary, since many functions lead to the same 
physical state. We may, for example, consider the case in which V'a-i anc ^ ^a-i are Slater 
determinants and Vyi-i = *&\-v ^ we then take tp(r) to be orthogonal to the orbitals in 
V'a-ij the antisymmetry requirement will turn the product into a Slater determinant and 
we obtain (f>i(r) = </>f (r)/\AA = ip(r). On the other hand, if tp(v) is not orthogonal to 
the occupied orbitals, the non-orthogonal components will be projected out as well. In 
the extreme case of <p(r) being a linear combination of the occupied states, 4>i(r) is zero. 
From these considerations we conclude that the difference between ip(r) and the associated 
physical state has no particular significance. The effect that antisymmetrization has on 
<f(r) is not pertinent to the one-body models considered in this paper since both already 
include antisymmetry, at least approximately. (The effect of antisymmetry on the shape of 
the spectroscopic amplitudes is distinct from the effect of antisymmetry contained in the N n 
of eq. (E3).) A more useful comparison would be between the spectroscopic amplitude and 
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(A/ii) 1 / 2 (^(r ,/ ), since this would measure the influence of the off-diagonal matrix elements of TV 



and thus test the validity of the single-particle model based on eq. [H)]). Neither comparison 
clarifies the role of the spectroscopic factors since, as previously noted, the spectroscopic 
factor is contained in both Xc and the spectroscopic amplitudes. 

At this point it is useful to return to the asymptotic normalization, A nc , of the spectro- 
scopic amplitude. For the low-energy 7 Be(p, 7) 8 B and 16 0(jo, 7) 17 F* reactions, for example, 
the capture occurs at large radii. Thus the asymptotic normalization of the spectroscopic 
amplitude is sufficient for describing the bound state in the physically relevant region. This 
is in line with the conclusions of refs. |8|, [|. As explained above, there are two different 
physical effects included in the spectroscopic amplitude and hence in its asymptotic normal- 
ization: one is related to the distortions of the (A — l)-body cluster due to the presence of an 
additional nucleon and is contained in the spectroscopic factor and the other is related to the 
single-particle properties of that extra particle and is described by the spatial dependence 
of the amplitude. Both are needed. 



V. CONCLUSIONS 

Spectroscopic amplitudes play a central role in the description of single-particle transfer 
reactions such as radiative nucleon capture. The amplitudes contain both single-particle 
and many-nucleon aspects of the nuclear many-body problem and can be, in principle, ob- 
tained from a fully microscopic model. We have presented two alternative approaches, each 
based on a set of coupled-channels equations, and have shown how one-body approximations 
can be derived in a systematic manner. We obtained two different single-particle models. 
In both cases, the single-particle wave function was found to be an approximation to the 
spectroscopic amplitude but normalized to one rather than to the spectroscopic factor. The 
quality of the one-body approximation depends on how well it describes the spatial depen- 
dence of the spectroscopic amplitude. However, even the best single-particle wave function 
will miss the spectroscopic factor. 

The first one-body approximation considered here results in A noninteracting particles in 
a single-particle potential. The simplification occurs since the equations for the spectroscopic 
amplitudes are no longer coupled; instead, we have A independent equations. Since we still 
have A particles and A orbitals, we can construct an antisymmetric A-body wave function 
by taking it to be a Slater determinant. In this model, it is redundant to explicitly enforce 
antisymmetry, e.g. through projection operators. The quality of this approach depends 
on how well the single-particle potential is chosen and how well the single-particle wave 
functions reproduce the shape of the spectroscopic amplitudes. 

The second one-body model presented here is usually derived within the generator- 
coordinate formalism. We obtained it by truncating a set of coupled integral equations. In 
this approach, antisymmetry is (approximately) imposed for each orbital separately through 
(approximate) projection operators. Unlike the first one-body approximation it can include 
many-particle correlations in the (A — l)-body subsystem. In principle, this approach is 
also simpler since we do not need to consider all the orbitals together. The price we pay is 
that the potential cannot be approximated as simply since it implicitly contains a projec- 
tion into antisymmetric states. There is also an explicit projection operator that must be 
approximated. The quality of this approach depends on the choice of the (A — l)-body wave 
functions and the importance of the channel coupling. 

It is important to realize that there is not one unique single-particle model. Different 
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one-body approximations to the nuclear many-body problem exist and can be derived in- 
dependently of each other. The resulting single-particle models may differ in subtle but 
crucial details and should therefore not be confused with each other. We have derived two 
such models and discussed their relation to each other and distinctions between them. In 
particular, the very different techniques for including antisymmetry should be noted. 

The spectroscopic factor, in its simplest form, reflects the partition probability of the 
A-body system into smaller clusters with allowance for antisymmetry effects. In the present 
paper it is cast in a complimentary light. It presents itself as a manifestation of the distortion 
of the (A — l)-body system due to the presence of the Ath particle. This distortion is 
both dynamical due to the interactions and kinematical due to antisymmetry. Since these 
are pure many-body effects they are, by definition, absent from one-body approximations. 
When center-of-mass corrections are included, the spectroscopic factors can be greater than 
one without violating the Pauli principle. Otherwise they must be less than or equal to one. 

A full calculation must include both the one-body and many-body effects. As the present 
paper emphasizes, this can be accomplished through the use of spectroscopic amplitudes. 
Obtaining these amplitudes, however, requires the full solution of the coupled equations 
presented here or a fully microscopic model. While this is still a distant goal, some recent 



fully microscopic approaches show promising results for low- mass systems [22 1 . Furthermore, 
models such as the shell model, the continuum shell model, or the cluster model include some 
many-nucleon correlations and provide reasonable approximations to the full problem. Since 
many-body effects are not contained in one-body models, cross sections calculated in this 
framework need to include the spectroscopic factor. For processes that are strongly peaked in 
the tail region, like the low-energy 7 Be(p, 7) 8 B and 16 0(jo, 7) 17 F* reactions, the one-body and 
many-body effects can be combined into a single parameter — the asymptotic normalization 
coefficient. 
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APPENDIX A: THE BASES 

We begin by defining a basis {^ A _ 1 } n =i,2,... f° r H^_ l5 the space of completely antisym- 
metric (A — l)-body wave functions ipA-i( r i, • • • , r A~i)- The basis states are orthonormal, 



/A-l 
\\ dr i ^X_ 1 (r 1 , ■ ■ • , r A -i)^A-i(n, ■■-, r A -i) = <W , 



(Al) 



and complete in H^_ 1 , that is: 

A—\ oo 



J] dr[ J2 n-iin, ■■■, r A -i)n*-i(r[, • • • , r' A _ x ) ^-xW, • • • , r' A _ ± ) 

i=l n=l 



if>A-i(ru--,r A -i) (A2) 
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holds for any ipA-i £ ^a-i- Specifically, for the sake of convenience, we choose \PX-i which 

are eigenstates of the (A — l)-body Hamiltonian H a-i = — Si=L 2m" + 2 Si 1^1 ^(l r * ~~ r j I)- 
To do so, we have to include both bound and scattering states. The superscript n labels the 
discrete as well as the continuous spectrum of H A _i. 

We now consider two different A-body spaces. The first, denoted by H^, is spanned by 

*7(n, • • • , r A ) = n-i(ri, • • • , ^^(r - r A ) , (A3) 

where r is a continuous parameter. The ^ n A are ortho normal with respect to both n and r: 

r A 

/ J] dr&Tiri, • • • , r A )* n A y (r l7 ■■■,r A ) = S nn ,S(r - r') , (A4) 

J i=l 

and the completeness condition for this basis is given by: 

J i=l \ra=l J / 

= il> A (r ir --,r A ) , (A5) 

where ip A G H^. The space H^ is a direct sum of the subspaces H^ and H^ which contain, 
respectively, totally antisymmetric and mixed-symmetry A-body states. The latter are an- 
tisymmetric in the first A-l coordinates and symmetric with respect to exchanges between 
the A-th nucleon and any other particle. For A— 2, the space H^l 2 is completely symmetric. 
An arbitrary wave function ip A G H^ can thus be written as the sum of an antisymmetric 
and a mixed-symmetric component: 

A ( A \ 

^(n, • • • , r A ) = —=i) A {r X i ■ ■ ■ , r A ) + ( 1 - —= \ i> A (ri, ■■■,r A ) (A6) 



= ^(r 1 ,---,r A )+K(ri,-'-,r A ), (A7) 

where ^ G H^ and ^^ G H^ 4 , and ^4 denotes an 'intercluster' antisymmetrization operator, 
which antisymmetrizes between the A-th coordinate and the remaining A-l coordinates. A 

is normalized by the condition A 2 = \J~A~A- Since (A /VA) (l — A/Va) = 0, the two 

subspaces are orthogonal to each other, that is: 

r A 

l[dr i if>?(r 1 ,---,r A )iftf(r 1 ,---,r A ) =0. (A8) 

"* i=l 

Furthermore, a symmetric operator O s cannot connect the two subspaces: 

r A 
J Hdniifin,- ■ ■ ,r A ) 6 s ^(n,- ■ ■ ,r A ) = (A9) 



i=l 



This includes the case where O s = H A , where H A is the A-ho&y Hamiltonian. 
The space H^ is spanned by 

A * n /(ri, ■■-,r A ) = A*Y = A[n-i(ri, • • • , r A ^)6(r - r A )} (A10) 
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or, equivalent ly, 

U*7) = a\r) \n~i) , (AH) 

where a^(r) creates a nucleoli at position r and we have used the convention ip A {r\, ' ■ ' > t a) — 
(7*1, ■ ■ ■ , r A \ip A ). The creation and annihilation operators a)(r) and a(r) obey the usual 
anticommutation relations, which ensure that the right-hand side of eq. (|A11|) is totally 
antisymmetric. The completeness condition for the basis {^^^ ,r } takes the form: 

/ A dr 'i ( T E / dr ^*A r (ri, • • • , r A ) A n r *(r[, • • • , r' A ) ) ^W, • • • , rJJ (A12) 

= iI)a(t x ,--,ta) , 



where ip A is a fully antisymmetric A-body wave function from H^. When ip A in eq. (|A12|) 
is replaced by a state from H^, the right-hand side of the equation vanishes. Thus V A = 
A~Y^=i$dr |_4\l/^' r )(_4\l/^' r | is a projection operator which projects states ipA £ EU onto 
their antisymmetric component ip A e H A . That V A = V A holds can be shown by using 
eq. ( |A14| ) below. 



The advantages of using totally antisymmetric basis states are obvious. The disadvan- 
tages of employing this basis lies in the fact that the states are no longer orthonormal. 
Instead, we have: 

r A 
\[dr iA ^Y\r 1) ---,r A ) A ^Y{r ir ..,r A )=M{n,r,n'y). (A13) 

"' i=l 

The norm operator, A/", and its kernel, J\f(n,r,n',r f ), have various interesting properties 
and are discussed in Appendix |B[ 

The basis states A ^ A r are not linearly independent, but are related to each other through 
the norm operator, J\ , as follows: 

1 °° r 

A n r (rx,---,r A ) = jJ2 drW(n,r,n',r') A ^ r '(r ir ..,r A ) . (A14) 

In fact, the {^^/^' r } basis is overcomplete and, at least in the case where the ^ A _i are Slater 
determinants, spans the space H^ A times. This accounts for the factor — in the first line 
of eq. ([l|), and in the completeness relation, eq. (|A12|) . 



An arbitrary antisymmetric A-body wave function i/i A (ri, • • • , t a) can now be expanded 
in one of the above bases. Using the set {\l/^' r }, we have: 

oo 

Mri, ■ ■ ■ , r A ) = -= J2 / dr *7(n, • • • , r A )Mr) , (A15) 

* n=l 

with expansion coefficients 

r A 

<f> n {r) = VI / H dr % *Y*(ri> ' " " > r A )Mn, ■■■,r A ). (A16) 

J i=\ 
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Alternatively, we can use the antisymmetric basis, {_4\l/^' r }, to write: 



Mri 



where 



1 °° f 

n=\ J 

A 

/ JJ dri A ^A r *( r ^ • ■ ■ ' r A )ipA{ri, ■ ■ ■ , r A ) ■ 



(A17) 



(A18) 



Equation ( A17 ) follows by applying the antisymmetrization operator to eq. ( A15|) and 
eq. (|A18|) can be derived from eq. (|A16|) by using the identity ip A (ri, • • • >t\a) = 
(A/yA) iPa{ti, " " " > r A ), which holds for totally antisymmetric A-body states, and the Her- 
mitean properties of A. Thus, the coefficients are the same in both expansions. 
From eq. ( |A15| ), it follows that is also possible to write the wave function if; A as 



Mri 



_. oo 
* n=l 



(p n (r A ) 



(A19) 



The coefficients <p n (f) are identical to those in the previous expansions. From eqs. ( All ) 
and (|A18|), one infers that they take the following form: 



4> n (r A ) = (^_!Kr)|^) . 



(A20) 



When ■0a(t , i ) • • • , r A ) denotes a bound state and the ^a-v which occur in the definitions 
of both \l/^' r and a® a r > are eigenstates of the (A — l)-body system, then the 4> n (r) are the 
spectroscopic amplitudes and the associated integrals S n = J dr \(f> n (r)\ are the spectroscopic 
factors. The A-dependent normalization factors in the above equations are included so that 
I lL=i dri \ipA(ri,-", r A ) \ = 1 holds, as well as 



oo „ 

J2 dr\Ur)\ 2 = A 

n=l J 



(A21) 



in accordance with the conventional normalization of the spectroscopic factors. This last 
equation follows by squaring eq. ( |A15| ), integrating over the coordinates, and using the 
completeness of the \l/^' r . 



APPENDIX B: THE NORM OPERATOR 

The norm operator for the antisymmetric basis {^^I / ^' r }, A/", and its kernel, J\f(n, r, n', r'), 
have many interesting properties. To start with, the kernel can be written in several equiv- 
alent forms: 



J\f(n, r,n, r') 



r A 

/ JJ dr t A^Y*( r h ■■■, rA)A^ A ' r '( r ^ ■■■,rA) 

J i=l 

VI / J] d Ti *T (n, • • ■ , r A ) A ^ r \ ri , ■ ■ ■ , r A ) 



(Bl) 
(B2) 



18 



(B3) 



= ^ ff[dr l ^Y*(r 1 ,---,r A )A^' r '(r 1 ,---,r A ) 

J i=l 

= Ajf[dr i n r \^---,r A )^=^ r \r u ...,T A ) (B4) 

= (*:U|a(r)oV)|*li) ■ (B5) 

We see that M{n, r, n', r') is not only the kernel of the norm operator for {^^/^' r }, eq. (|B1|) , 
but is also proportional to the overlap of an element from {yi^/^' r } with an element from the 
non-antisymmetrized basis {\l/^' r }, eq. (p2|). Furthermore, N{n,r,n',r')/v~A is the matrix 
element of the 'intercluster' antisymmetrization operator, A, in the basis {\l/^' r }, eq. (|B3|) . 
or - equivalently - J\f(n, r, n', r')/A is the matrix element of the projection operator A/yA 
in that same basis, eq. ( |B4| ). Finally, we can write J\f(n,r,n',r') as the matrix element of 
a^a^ir') in the (y4-l)-body basis, eq. ( |B5| ), where a^ir) and a{r) create and annihilate, 
respectively, a nucleon at position r. 

Since J\f(n, r, n', r') is proportional to a projection operator, it has no inverse. However, 
its square root - in the sense of a matrix operation - can be given. It is simply the matrix 
element of the antisymmetrization operator A: 

i r A 

y/tf(n, r, n', ?) = J J[ dr, tf 7*(n, • • • , r A )A^Y (n, • • • , r A ) ■ (B6) 



To see the projection operator nature of N more directly, one can multiply eq. ( A15 ) by 
A^ n A r (fi, • • • , T A ) and integrate over the coordinates. If ip A is completely antisymmetric, 
we obtain the following equation for the expansion coefficients: 

oo „ 

Mr) = jJ2 ^r' JV>,r,n>')<MO . (B7) 

n'=l J 

For a mixed-symmetric state ip A , on the other hand, we find: 

1 °° f 

o = -j y. / dr ' u ^ r > n '> r ')^' ( r ') > ( B8 ) 

i.e. J\/A, when acting on a set of expansion coefficients <p n { r )i behaves like a projection 
operator: It returns the coefficients <fi n { r ) of an antisymmetric state and yields zero when 
the 4> n { r ) correspond to a mixed-symmetric state. It follows immediately that in a restricted 
space of coefficients which originate from a completely antisymmetric wave function, J\/A 
becomes the identity matrix. In this restricted subspace, the functions A ^> n A ' r /VA act in 
many ways as if they were orthonormal. 

If we expand an excited state of the A-body system, ty A (ri, • • ■ , r A ), as in eq. (|l|), we 
obtain expansion coefficients (p^ ,k (r) = (\l/^_ 1 |a(r)|\l/^). For k — 1, ^ A =1 describes the 
ground state of the A- nucleon system and the 4>n' k=1 { r ) reduce to the usual spectroscopic 
amplitudes. For a fixed n, on the other hand, the 4>n ,k (r) correspond to particle states built 
on |\E f ^_ 1 ). Similarly, one can expand an (A-l)-body state, N^X-n in terms of (v4-2)-body 
basis states, ^ A _ 2 , and obtains expansion coefficients 4>^ 1,n {r) = (\l/™_ 2 |a(r)|\I/^_ 1 ). With 
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respect to the (A-l)-body system, the (A-2)-body functions represent hole states and the 
A-body functions are particle states. The expansion coefficients 4>^~ 1 ' n {r) can be used to 
rewrite the equations of motions (eq. (|)): 

(E A - El_ x )4> n {r) = -^Mr) 

CO CO „ 

+ E E / ^^V^'-VMk' - r\)Mr) ■ (B9) 

n'=l m=\ •* 

The kernel of the norm operator can also be written in terms of the particle or hole 
states. For the particle states, we insert a complete set of intermediate A-body states, 
^(r*!, • • • , r A ) in eq. ( |B5| ) to obtain: 

OO 

AT(n, r, n', r') = £ ^ k (r)^ M (r') . (BIO) 

fe=i 



On the other hand, using eq. ( |A3| ) , the kernel of the norm operator can be expressed in 



terms of the (A-l)-body wave functions: 
J\f(n,r,n',r') = 5 nn >5(r - r') - 



„A-2 

(A - 1) / H dr^ n /-i(rx, • • • , r A - 2 , 0*li(n, • • • , r A _ 2 , r) . (Bll) 



This is not diagonal, in either n or r, even for the simplest systems. From the last expression 
we derive the hole-state form: 

J\f(n, r, n\ r') = 5 nn >5(r - r') - 

OO 

E^~ 1,n v)^- 1,n '(»o- ( B12 ) 



771=1 



Combining eqs. ( |B10| ) and ( |B12| ), we obtain a completeness relation for the spectroscopic 
amplitudes corresponding to the set of particle and hole states: 

CO CO 

5nn,5{r-r') = £ ^^(r 1 )^ 1 ^ (r) + $>^ fc (r)<^V) • (B13) 

m=l fc=l 

The sum over m runs over all states of the (A — 2)-body system [the hole states of the (A — l)- 
body system] while the sum over k runs over the states of the A-body system [the particle 
states of the (A — l)-body system]. The spectroscopic amplitudes for the particle states 
are not complete by themselves since they lack the contributions that are Pauli blocked, 
namely those contributions corresponding to hole states. Contrary to the impression that 
this equation may give, the 0's are neither orthogonal nor normalized to one. 

To illustrate the formalism, we consider a two-particle system. In this case, the 
^A-i( r i5 • • • 5 r A-i) = ^i(r) are one-body wave functions and the kernel of the norm oper- 
ator is given by: 

M{n, r, n', r') = 6 nn ,6(r - r') - ^*(r')^i(r) . (B14) 
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The terms diagonal in n are projection operators onto states orthogonal to \l/"(r). This is 
also true for larger particle numbers if the \I/^_ 1 (r 1 , . . . , Va-i) are single Slater determinants 
[In that case Af(n, r, n, r') = S(r — r') — Y^ s =i ( t ) T( r ') ( f ), s ( r )> where the sum is over occupied 
single-particle orbitals <f>™(r')]. While the matrix elements diagonal in n are projection 
operators by themselves, one needs to divide the full norm operator by A in order to obtain 
a projection operator. 

APPENDIX C: A BOUND FOR THE SPECTROSCOPIC FACTOR 

In this appendix, we show that S n can be written as the product of two factors, which 
express antisymmetry and dynamic distortion effects, respectively. When recoil and center- 
of-mass corrections are neglected, both factors have to be less than or equal to one, yielding 
an upper limit of one for the spectroscopic factor. 

We start by defining a normalized spectroscopic amplitude 0„(r) = <f> n (r) / yfS^ and 
express y/S^ as follows: 



<S n = J dr<f>* n (r)<f> n (r) = (A[<f> n n-i\M • (CI) 

Next, we introduce a projection operator 

_ \A[4> n n-i])(A^n-i]\ (r9 , 

^n = 77 , O^J 

where N n = (A[<f> n ty A _ 1 ]\A[<f> n ' i &' A _ 1 ]). An arbitrary state \ipjy) can then be broken into two 
orthogonal parts: 

Wa) = Pu\^a) + (1 - PuMa) = WY) + \^ n ) , (C3) 

where \if> A ) and \ip A ' n ) are the components of \i/)a) which are parallel and orthogonal, 
respectively, to the state ^[^n^l-i])- We can then write the spectroscopic factor as: 

S n = N n {i>YW)f) = N n ((^a) - ^i' B |^i ,B >) • (C4) 

The expression in brackets is less than or equal to one since (i^a\^a) — 1 an d both {VaI^a) 
and (if) ' A n \if ' A ' n ) are positive semi definite. When the (A — l)-body system is completely 
described by the wave function ^ / ^„ 1 , i.e. when there are no distortions due to the potential 

of the A-th. nucleon, (if> A n \ip A ' n ) vanishes and (ip A \if> A ) = 1. Thus (if) A \if) A ), and therefore 
S n , provides a measure of the dynamic distortions induced by the presence of the extra 
particle. 

Since the factor N n can be expressed as: 



oo / „ x 2 



ra=l 



N » = 1 ~ E / dr^t^MMr) , (C5) 



where the sum is explicitly non-negative and less than or equal to one, it is also restricted, 
N n < 1. N n carries the effect of the antisymmetrization and equals one only when <f> n (r) is 
orthogonal to <f ) ^ 1 ' n ( r ) f° r a U m - 
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From the above considerations it follows that S n < 1. If antisymmetry was neglected, the 
spectroscopic factor could be as large as A, since the sum rule given in eq. ( |A21| ) would be 
the only restriction on S n . When center-of-mass corrections are incorporated, eq. flU4|) still 
holds, but eq. (|C5[) has to be modified and N n can become larger than one. The influence 



of center-of-mass corrections is discussed in the next Appendix. 

APPENDIX D: CENTER-OF-MASS CORRECTIONS AND INTRINSIC 
SPECTROSCOPIC AMPLITUDES 

When dealing with the center-of-mass problem it is useful to introduce the Jacobi coor- 
dinates Pj = Rj — rj+i, where Rj is the center-of-mass coordinate of the j-body system 
defined by particles 1 through j. Taking into account the center-of-mass motion, the A and 
(A — l)-body wave functions are written as: 

exp[ik A - R A ] , 
4> A (r 1 ,---,r A ) = l —j= ^i(pi,"-,PA-i) (Dl) 

and 



V 



«^-'(n, • ■ • , rx-0 = °Mik A j_R A ^ ?{pu . . . pA2) (D2) 

respectively. Here k A and k A _\ are the center-of-mass momenta of the A and (A — l)-body 
systems, respectively. We have used box normalization with volume V. The spectroscopic 
amplitude is written as: 

A 

<Pn,k A -Ar) = VA ]Jdr i 6(r-r A )'$ n > kA -i*(r 1 ,---,rA-i)Mrir--,rA) (D3) 

expfir ■ (k A ~ k A -i)] rj f rr , ex P [iPA-i ' (^r fc A - k A -i)] 

= vv ^J l\ dpl vv 

xV?(p!, • ■ • , p A - 2 )^ A (pi, • • • , Pa-i) (D4) 

eM^r-(k A -k A ^] 2l fA-l kA _ k \ (D5) 



Vv n V a 



This equation is unexpected and requires some comments. Formally it is correct: the spatial 
dependence of the spectroscopic amplitude is given by a plane wave and the spectroscopic 

2 

factor is <^ (^j^-k A — fc^-i) • Since 4> n {k) is on the order of 1/vV, it is small and the 

condition that the spectroscopic factor must be less then or equal to one is easily satisfied. 
The plane wave behavior of the spectroscopic amplitude arises from translational invariance. 
The combination (^j^-k A — fe^-i) is Galilean invariant. By taking both the (A — l)-body 
and the A-body systems to be in states of good momentum we have forced the Ath particle 
to also be in a state of good momentum; 0^ (k) is then the probability amplitude for finding 
the Ath particle with relative momentum k when the (A — l)-body system is in state n. Its 
Fourier transform, which we identify as the intrinsic spectroscopic amplitude, is given by 
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(compare eq. ( |A16| )): 



a 'p) = I dk A J Mlp ^ A - l] ~^ n {k A _ l )= I dk A ^ n , kA _ x {p) 



k A =0 



r A-l 

VA Y[dpi S(p - p A -i)VT(pi, ■■■, Pa-2)iI>a(Pi, ■iPA- 
J 1=1 



(D6) 
(D7) 



In analogy with 0^ (k), 0^(p) is the probability amplitude for finding the Ath particle at 
the distance p from the center-of-mass of the (A — l)-body system when that system is in 
the state n. The intrinsic A-body wave function can be written in terms of the intrinsic 
spectroscopic amplitudes as (compare eq. ([[])): 

^(pi, • • • , PA-l) = JjJ2 0»(PA-1)*/(P1, • • • , PA-2). (D8) 

* n=l 

As we show in the next paragraph, the intrinsic spectroscopic amplitude is also the quantity 
that is needed to calculate physical observables. 

We now write the transition matrix element, eq. (|T8|), in terms of the intrinsic spectro- 
scopic amplitude as follows: 

A oo „ 

<Vu|^expHfe.r,]|V£> = E J dr( t>nM-M) ™pH* • H< feA _» (D9) 



8=1 n.,&A 



J2 d P <f>n(P) ex PH^ " P(A - l)/A]cj>«\p) 

n=l ^ 

x<5(A; x + A;-A; A ) , (D10) 



where the transition operator has been taken to be a plane wave as is appropriate for 
radiative capture and the spin and isospin dependencies have been suppressed for simplicity. 
The functions 0^ 7 (p) and 4>n{p) are intrinsic spectroscopic amplitudes for the scattering 
and bound states, respectively. The delta function ensures overall momentum conservation. 
The (A — I) /A factor in the exponential takes care of the laboratory to center-of-mass 
transformation. 

The equations of motion for the intrinsic spectroscopic amplitudes are easily derived by 



substituting eq. (fD5|) in eq. (p9| ). This gives: 

vv2 ^ 



(£% - El-i)<f£(p) = -^(P) + EE dp'p nn ,{p')V{\p'-p\)4>l{p), (Dll) 
where \i is the reduced mass and 

Pnw(p) = {j^2J ^(PiA - l)/(A - 2))tf? \p(A - I) I {A - 2)) (D12) 

is the transition density for the (A — l)-body system. The (A — I) /(A — 2) factors originate 
in the conversion from the Ra-2 — r A-i coordinate to the Ra-i — Va-\ coordinate. The 
diagonal transition density is the usual density and is normalized to A — 1. 



2:3 



The remaining quantity to consider is the norm operator, Af '. This is most easily done 
starting with eq. (|B12| ). We obtain the following expression for the kernel J\f(n, r, n', r'): 



Ar(n,r,n',r') = 5 nn ,8(r - r') - ( it ' 



,2 \3 



oo 
X 

in 



\A(A-2] 
V 3r ( {A - 1)2 (r> 



(A -I) 



In contrast to the situation where the center-of-mass corrections are neglected, the am- 
plitudes in the sum given here depend on both coordinates. Consequently, the intrinsic 
spectroscopic factors S^ = J dp \<j>n(p)\ no longer have to be less than one. This is il- 



lustrated in ref. |L5] for the harmonic oscillator model. The completeness relation for the 



particle and hole states, eq. ( P13|) , is also modified, since this last equation must be used 



instead of eq. ( |B12|) . The spectroscopic amplitudes corresponding to the particle states are 
just replaced by their intrinsic counterparts. 
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